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J> ' Abstract. We present new analytical results concerning the spectral distributions 

■ for (2 x 2) random real symmetric matrices which generalise the Wigner surmise. 
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1. Introduction 

The level statistics of a quantum system represents the most significant, although not 
the only M, signature of quantum chaos. The Poisson and Wigner distributions of 
k*" ■ dimensionless nearest-neighbour spacing, s, 

p P (s) = exp (s) , (1) 

Pw(s) = yexp(-7rs 2 /4) , (2) 

are known in quantum chaos theory as two universalities that correspond to two extreme 
cases of classical dynamics, namely purely regular and completely chaotic (see , e.g., 
0). The majority of many-body systems such as nuclei, molecules, atoms or solids (see 
|1|, |^, U and references therein) have been found to be chaotic although for such complex 
systems no classical limit can be constructed. 

As was first recognised by Wigner 0, the nearest-neighbour spacing distribution 
(NNSD) (Q) well corresponds to the eigenvalue distributions of random matrices, and 
this explains the importance of the Random Matrix Theory || for studying statistical 
properties of many-body systems. Particular attention was paid to Gaussian ensembles. 
In fact, assuming that (i) the elements of the Hamiltonian matrix are independent real 
variables and (ii) the matrix distribution is invariant under an orthogonal transformation 
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of the basis states (see, e.g., chapter 3 in 0), one finds that the matrix elements 
are independent Gaussian variables with zero mean and with variance satisfying the 
conditions afj = (1 + 5ij)a 2 . Imposing particular symmetries on the Hamiltonian, one 
gets J7| Gaussian Orthogonal, Unitary or Symplectic Ensembles (GOE, GUE or GSE, 
respectively), which are widely and successfully applied in many fields of physics (see 
e.g. reviews |, |, |§). 

Up to now, excited atomic nuclei are considered to be the best examples of chaotic 
quantum systems 0. Starting from slow-neutron scattering experiments, which were 
first described in terms of random matrices ||, a lot of nuclear structure data has 
been analysed in the context of chaos. However, it was repeatedly noticed that 
the experimental data do not exactly match the distribution (|2]) and exhibit slight 
deviations [§, ||. These deviations are thought to be caused by the fact that the real 
system is not purely chaotic, but can be the quantum analog of a classical system that 
is transitional between chaotic and integrable. In this context, a few phenomenological 
formulae were proposed and analysed. The most famous are the Brody |§ distribution, 

p w (s) = (u + l)as» exp (-c^ +1 ) , a = [r )] + , (3) 

which was shown to match better the experimental data || on both high- and low-energy 
nuclear spectra, and the Berry-Robnik distribution [|TD| 

Pbr(s) = e( ,? - 1 ) s {(l-g) 2 erfc( v ^ : g S /2) + [2g(l-g) + (7r/2)g 3 S ]e-^ 2s2 } .(4) 

Although the Brody distribution takes the form of Poisson for u = and Wigner for 
L) — 1, it has the unrealistic property that its derivative goes to infinity at e = |TT|] . 
Moreover, the parameter u has no clear physical meaning. On the other hand, the 
distribution @ does not give a level-repulsion for non-integrable systems. Caurier et 
al [11 considered a model which allowed to simulate the transition from integrability 
to chaos and succeeded to derive the asymptotic limits for small and large neighbour 
spacings in the near-integrable limit. 

The idea to derive a third universality class corresponding to intermediate statistics 
has been actively pursued in recent years. In particular, in Refs. [|l|, 13[], the distribution 
of n nearest neighbours, 

(ft -L 1 

P m (n, s) = Yj^TTF * ^ 6XP + 1)S] ' (5) 
was shown to be relevant to a certain class of exactly solvable models with nearest and 
next-to-nearest neighbour interactions, generalizing the results obtained earlier |14|, pl| 
on pseudointegrable billiards and the short-range Dyson models. Remarkably, the 
NNSD given by (||) with n — 1 exhibits a level repulsion ~ s 13 and falls to zero at 
s as exp [— (/3 + l)s]. For (3 — 1 it is referred to as the semi- Poisson distribution jj"4l |I5(|. 

In this letter we argue that one might search for an explanation of the discrepancy 
between data and random matrix theory simply by generalising the random matrix 
ensemble used. In fact, the Gaussian ensemble is defined by the two assumptions 
(i) and (ii) mentioned above, and their applicability should be carefully checked for 
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a given physical system. As far as nuclear physics is concerned, the deviations of 
the experimental data on slow-neutron or (p,p') resonances from a random matrix 
description could simply arise from the non-invariance of the random matrix ensemble 
under an orthogonal transformation of a basis, i.e. the assumption (ii) is violated. 

In addition, it is well known that realistic interactions in many-body nuclear, 
molecular or atomic systems are predominantly of one- and two-body nature, implying 
that the distribution of the matrix is not only not invariant under an orthogonal 
(unitary) transformation of the basis, but also that the elements of the Hamiltonian 
matrix are not independent (if the number of particles is more than two), i.e. both 
assumptions for a Gaussian ensemble do not hold any longer. In this context, French and 
Wong |l6j and Bohigas and Flores [jl7[] independently introduced the Two-Body Random 



Ensemble (TBRE), which is characterised by a Gaussian level density distribution, 
rather than a semi-circle provided by a GOE ||16|| . The level distributions of experimental 
nuclear spectra did confirm this result. The NNSD relevant for a TBRE was found 
numerically to be fairly represented by the Wigner surmise although recently it has 
been pointed out |19| that NNSD given by (|5|) with n = 1 and a certain real value of j3 
fits better the shell-model spectrum obtained with the sd-interaction of Wildenthal |20[ . 



Given the importance and actuality of these investigations, we present in this letter 
some analytical results concerning the properties of (2x2) random symmetric matrices 
for which the assumption (ii) mentioned above is not satisfied. First, we derive the 
Hamiltonian distribution as a function of its eigenvalues and we calculate the NNSD 
which generalises the well-known Wigner surmise ||. We show that the model allows 
to describe the transition from purely chaotic to asymptotically nearly integrable limits, 
being different from the intermediate statistics mentioned above. Then, for a particular 
case, we give the analytical expressions for the lowest moments of this distribution. 
Finally, we propose a method to derive the moments of the eigenvalue distribution 
without knowledge of an explicit expression for the distribution. 

2. Generalised Wigner surmise 

Let us consider a (2x2) real symmetric matrix 

H=(" U f 1 ), (6) 

V n 21 -"22 / 

whose elements are independent Gaussian variables with zero mean and variance afj, 
and if 12 = H^i- The probability density of the matrix H is then given by 



p(H) = ) exp 



tt2 tt2 tt2 

11 _|_ 12 _j_ 22 



2al 2a?o 2a? 



(7) 



'11 ^12 ^22, 

Each matrix H can be diagonalised in an orthogonal basis and therefore H = O t DO, 
with 

/ cos 9 — sin 9 \ 
\ sin 9 cos 9 / ' 
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and 



D 









Ep 



Similar to the case of GOE pi ], we find that in the general case 



H22 



E a cos 2 9 + Ep sin 2 9 
(E a — Ep) cos 9 sin 9 
E a sin 2 9 + Ep cos 2 9 



We deduce that the probability density expressed in terms of the eigenvalues and the 
angle 9 is 



p{E a ,Ep,9) 



En, — Eh 



(2vr) 3 / 2 V^i 



rr 2 rr 2 
°22°\2 



■ exp 



E a T? - (E a - Ep) (a 2 u cos 2 9 



a 22 sin 2 9 



2 \ 



2a 2 1 a 2 2 £ 2 



cxp 



-~(E a - Ep) 



'cos 2 (29) sin 2 (29)' 



Aa 2 



.1.2 



(8) 



where E 2 = o\ x + o\ 2 and E a — Ep > 0. 

The nearest-neighbour spacing distribution for the variable e = E a — Ep 
with the mean spacing D = J ep (e) de, is given by the following integral 



Ds, 



p(e) 

from which we obtain 

p(e) 



7I-/2 

-7t/2 



d9 / dE n 



2^/S 2 ( T 2 2 



: CXp 



dEp p (E a , Ep, 9)S(e-E a + Ep) 
e 2 (S 2 + Aa 2 2 ) 



les 2 ^ 



Aa 2 2 ) 



(9) 



(10) 



12 



where Iq is a modified Bessel function of the first kind. 

The expression ( |TU| ) looks like a Rayleigh-Rice distribution, well known in signal 
theory p2 |, except for the argument of Iq, which is not linear as in the usual Rayleigh- 
Rice distribution but quadratic. This is why we will refer to p (e) as to a quadratic 
Rayleigh-Rice distribution. 

Let us consider a particular case when the diagonal matrix elements have the same 
variance o~\ x = a 22 , which is x times larger than the variance, a\ 2 = a 2 , of the non- 
diagonal matrix elements, i.e. x — cr ii/ cr i2- Then the eigenvalue distribution (H) reduces 
to 



p x (E a ,E / 



E n — Er~ 



7V 



exp 



E 2 a + E 2 + l(E a -Ep) 2 ( X -2)sm 2 (29) 



(2no- 2 fl 2 x 

while for the nearest-neighbour spacing we get 

(x + 2)e 2 



'ID 



Px( £ ) 



2 X 2a- 



cxp 



(x-2) g s 
16x<r 2 



(12) 



lQxo~ 2 J 

For x — 2, expression (|TJ) reduces to the Wigner surmise. The distributions p x (e) are 
plotted in figure 1 for x = 1, X = 2 and x = 5- 
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For small values of e, the distribution ( |1~2| ) goes linearly to zero, 

1 



Px( £ ) K /tt^o 2 g , (13) 



while for large e, 




^2 

(14) 

X = 2 



The asymptotic behaviour has a functional dependence on e similar to that of the Wigner 
surmise, i.e. the NNSD (|TT| ) goes linearly to zero for e — > and it falls down according 
to exp (— e 2 ) for e — ► oo. However, as seen from (|l3|)-(|i~4T), the natural dependence on 
X provides a certain scaling. 

To calculate various statistical characteristics, it is often required to know certain 
moments of the distribution. Thus, we have derived analytical expressions for some 
moments of the distribution (|1^), and the lowest are given in table [I]. 

From (|T2|) and the expression of the mean spacing, D (cf table [I]), we have derived 
the distribution of the dimensionless nearest-neighbour spacing, s = e/D : 

q x (s) = Dp x (sD) . (15) 

For x = 2, one finds the Wigner distribution (||). It can be shown that for \ an d 
x' = 4/x, the two functions q x and q x < are equal. 

The distributions q x for x = 4, x = 30, x — 500 and x — 1000 are plotted in 
figure 2(a-d) and compared to the Poisson (dashed), Wigner (dotted) and semi-Poisson 
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Table 1. The moments up to n = 5 of the quadratic Ray leigh- Rice distribution (Q) 
for er 2 = 1. For n = 1, one obtains the mean spacing D. i* 1 is the hypergeometric 
function | j23| . 



M„= f °°£"p x (e)de 



values for x = 2 for x = 5 



n = 1 
n=l 4xv / 2^(x + 2)- 3 / 2 J F 



2 (2 + x) 



= 3 96x 2 v / 2^(x + 2)~ 5 / 2 F 



5 1 i (X=l)' 
4' 4' x ' U+2 y 



71 = 4 



4 (12 + 4x + 3 X 2 



5 3840x 3 V2^(x + 2)- 7 / 2 J F 



Z 2 1 (x-lY 

4' 4' \^ x +2 J 



1 1 

\/27r 1.3 V2tt 

8 14 

12v^tt 28.7v / 2tt 

128 428 

240^ 1139. 2V2T 



(dashed-dotted) distributions. Note that the semi-Poisson distribution can be fairly 
approximated for x — 30 although its asymptotic behaviour for large s is different. As 
follows from figure 2, the parameter x allows to describe a transition of a quantum 
system from a completely chaotic limit (x = 2) to a nearly integrable one (x — > 00 or 
X — *■ 0). However, the integrability characterised by the Poisson distribution ([!]) can 



never be reached (compare with the model of Caurier et al |TT[| ). 

Indeed, integrating p x in (|Tl~|) over Ep from -00 to i? Q and over E a from —00 to 
+00, we obtain the angular distribution 

r- (6) = =r . (16) 

vrV2 [i+ i( x _ 2 )sin 2 (2^)] 1 ' 

This distribution is represented in figure 3 for different values of x- F° r X = 2 it is 
an exactly uniform distribution, which means that there is no privileged basis (the 
orthogonal invariance holds). For high values of x-, the initial basis is nearly the 
eigenbasis (the diagonal elements are much larger than the non-diagonal ones), thus 
r x takes its maximum absolute values for 6 = and 9 = n/2, whereas for small values 
of X-, the eigenstates are more likely obtained after a rotation of tt/A of the initial basis 
and r x is maximum for 9 = 71 /A. For x and x' = 4/x the two curves are in quadrature. 

We can re-express p x (E a , Ep) as a function of e and S = E a + Ep. Then p x can be 
factorised into a function depending on e times a function depending on S, i.e. these 
variables are independent. Moreover since S is the trace of the matrix it is a Gaussian 
variable with zero mean and all its odd moments are zero. From the independence of e 
and S we deduce that the moments of the eigenvalues fulfill 

(K) = (-i) n (^ n ) (17) 

(E n a ) = ^J2{2 P )(S 2p }(e n - 2p } (18) 



From (0) and (|I8"), we can derive the moments of the eigenvalues whose distributions 
are difficult to compute. Note that we deduce from fllTD that the highest and the lowest 
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Figure 2. Nearest neighbour spacing distributions q x (solid line) for \ = 4 (a), \ = 30 
(b), x — 500 (c) and \ = 1000 (d). The Poisson distribution is plotted in dashed line, 
the Wigner distribution in dotted line and the semi-Poisson in dashed-dotted line. 

eigenvalues have opposite mean values and the same variance. 
3. Conclusion 

The study of the statistical properties of spectra of realistic Hamiltonians requires the 
consideration of a random matrix ensemble whose elements are not independent or whose 
distribution is not invariant under orthogonal transformation of a chosen basis. In this 
letter we have concentrated on the properties of (2x2) real symmetric matrices whose 
elements are independent Gaussian variables with zero means but do not belong to the 
GOE. We have derived the distribution of eigenvalues for such a matrix, the NNSD 
which generalises the Wigner surmise and we have calculated some important moments. 
The asymptotic properties of the distribution obtained are functionally identical to 
those of the ordinary Wigner surmise. For finite \i the model considered here allows to 
describe the transition from chaos to near integrability (the exact integrable limit is never 
realized). Thus it represents a chaotic system although with a degree of disorder less 
important than in the Wigner surmise. The derivation of similar analytical expressions 
for matrices of larger dimensions is technically difficult. However, we believe that the 
present results already justify the use of NNND of type (|12"D to fit the data as an 
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6 (rod) 



Figure 3. Angular distributions for \ = 1 (dotted line), \ — 2 (dashed line) and 
X = 5 (solid line). 

alternative to the Brody distribution. We also think that these results can be extended 
to hermitian matrices. 
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